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1. Introduction 

The present paper will revisit the celebrated Berezinsky-Kosterlitz-Thouless (BKT) 
phase transition |^ aiming to address the quantum turbulence at two dimensions. For 
our purpose, we'll utilize the path-integaral formulation of the Abelian-duality at two 
dimensions, which is developed in . 

In 1973 two articles appeared at the same time. S. Coleman's paper shows that: 
the continuous symmetry of two-dimensional quantum field theory is not spontaneously 
broken. This is not difficult to prove: If a continuous symmetry spontaneously broken, 
there must be massless Nambu-Goldstone boson states in the physical spectrum, which 
is contradictory to that, in two dimensional quantum field theory, there are no any phys- 
ical local operator can create massless states. The latter is because, from the spectrum 
representation of two point correlation function, you can know that any two-dimensional 
massless states will contribute a non-positive definite item — Clog(x^), thereby undermin- 
ing the positivity of the two point function of physical local operator (so any such local 
operator must be un-physical) . 

In fact, in statistical physics, even a few years ago ^ Mermin, Wagner and Hohen- 
berg have discussed the corresponding conclusions of Coleman's theorem. In their discus- 
sions, thermodynamical fluctuation will replace the location of the quantum fluctuations, 
destroy any effects of continuous symmetry breaking. 

In addition, the well-known Landau theory of phase transitions tells people: continu- 
ous phase transition and spontaneous symmetry breaking is always closely linked. There- 
fore, according to people's understanding at that time, Coleman's theorem roughly means 
that for a two-dimensional system only with continuous symmetry, phase transition can 
not occur! 

In 1973, another interesting article study such a system appeared. In this classic 
article, Kosterlitz and DJ Thouless studied a SO{2) Heisenberg spin model on lattice 
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plane. They found a continuous phase transition! They found that this system has a 
critical temperature, when system temperature is lower than it there is a long range order, 
when higher than the critical temperature only shot range order can exist! (This system 
has been studied by Soviet physicist Berezinsky in 1970, but the article appeared in the 
Soviet Union's journals.) 

This seems contradict to S Coleman's theorem, what's wrong? who is right? 

Later, it was clear, S Coleman and BKT both right. Two-dimensional continuous 
symmetry is not spontaneously broken, while in the 5*0 (2) model, there is indeed a long 
range order to a short range order continuous phase transition. The key point is that 
the continuous phase transition is not necessarily accompanied by spontaneous symmetry 
breaking! 

In fact, the SO{2) model on two-dimensional plane has stable vertex excitations (with 
topological quantum numbers). At low temperatures vertex and anti- vertex pair up, those 
have little effect to the thermodynamic behavior of the system. One need only consider 
the phase fluctuations, the long range behavior of the system consists of a free-field theory 
description, between two phases (or between two SO{2) spins), there is a long-range cor- 
relation. However, when the temperature goes up to a critical point, inverse temperature 
(3 = 1/T times the energy required to unlock the corresponding vertex pair is lower than 
Entropy. Therefore, vertex and anti-vertex does not tend to twinning, the un-twinning of 
vertex anti-vertex pair causes the phase transition. Above the critical temperature, the 
gas of vertex and anti-vertex produced a Debye screening effect, undermine the long range 
order. At the critical point, the fluctuations of vertex is marginal, the system will flow to 
a non-trivial fixed point, associated with algebraic long-range order. 

2. KT Phase Transition and Quantum Turbulence 

2.1. Field Theory Discussion 

Since we only interest long-distance behavior of the system, we can develop a contin- 
uous field theory description to the BKT phase transition, to achieve this, we'll utilize the 
Abelian-duality for Abelian field theory at two dimensions, which is develop in 0] . 

Firstly, we rewrite the 5*0(2) spin field J(x) as J = exp(zA), where A[x) is a contin- 
uous scalar field. We ignore the vertexes at first, then the system partition function can 
be give as an appropriate path integral of A field, the weight of this path integral is. 
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W[A] = exp 
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(2.1) 



where T is temperature, dA is the differential of A, and | |c/^| P = {dA) A *{dA), * is Hodge 
star operation at two dimensions. 

Now, we will add the vertex contribution. When there is a vertex singularity at point 
p, we'll take an infinite small circle C around p, we cut out the inner of C - thus the vertex 
singularity is cut out - and replace the path integration of the field configuration over it 
by a small constant g, which marks the contribution of the inner of a vertex. Then we 
take the configuration of A field along C (which is determined by the vertex) as boundary 
condition of the field configuration at the outer. Finally, let us perform the path integral 
at the outer of C with vertex boundary condition. The total contribution of a vertex is 
the multiplication of this outer path integration and factor g . 

In fact, from the definition of a 1- vertex (which carries topological quantum number 
1), the configuration of A along C can be simply characterized as: The value of A will 
increase 2n after walking around C counterclockwise. If what we are considering is anti- 
vertex, then the value of should be reduced 27r, not increase. In another word, the value 
of A should increase — 27r for anti-vertex case. 

For the vertex with topological quantum number 2, 2-vertex, we can image that it 
is composed of two tightly closed 1-vertexes (to form a single 2-vertex). However, be- 
cause vertex - vertex is always mutually repulsion, so the two composite 1-vertex always 
tend to separate from each other, the separate configuration have lower energy and thus 
higher integration weight. Therefore, we'll firstly ignore all higher vertexes with topological 
quantum numbers greater than one, but only consider the contributions of well separated 
1-vertexes, and 1-anti- vertexes (-1-vertexes), they form a kind of dilute gas. 

In principle, for each of vertex dilute gas field configuration we can calculated the cor- 
responding path integral, the calculation is simply a generalization of the above calculation 
to single 1-vertex. We will add all the contributions, including the field configuration of 
no vertex, all add up, end up with the partition function of the system: 



the dilute gas of m 1-vertex and n -1-vertex, the position coordinates of the vertexes and 




(2.2) 



Where J[T>A]m,n stands for the path integration for the A field configuration with 
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anti-vertexes are also integrated. And the factor comes from the permutation of 

identical particles, factor g'^+'^ comes from the inner configurations of m vertexes and n 
anti-vertexes. 

So far, it seems all is well, the left work is to carry out the path integration with 
1-vertex boundary condition. 

2.2. Abelian Duality 

We begin from the configuration without vertex. In short, we imitate the Fourier 
transformation to transform the path integral of A field to path integral of some B field, 
by noticing that the path integration J[PA]1^[A] is equivalent to 



Z = [ [DAVB] exp — / {Fb ^ dA) 
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(2.3) 



where Fb = dB. One can see that the exponential factor is automatically trivial by 
integrating by part. Thus path integral over B only bring us a trivial infinite constant, 
which can regularized easily. 

Now we exchange the order of A, B functional integral. Note that we are considering 
is a topological trivial A field (without vertexes), the A functional integral and functional 
integral of Fj\ = dA is basically the same. Therefore, we first perform the integration of 
dA, easily calculated, the result is 



/ [PS] exp 
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(2.4) 



In other words, with out vertexes the partition function J[I'A]W[A] is completely equiv- 
alent to Z' (|2.4|) , in terms of functional integration over field B. 

Now, if you want to calculate, say /[I^A](W^[A](iA). The same transformation as 
above will tells us that it equivalent to calculate 



M](iT*dS)exp 



T 

47r 



\dB\ 



(2.5) 



With these preparations, we now can perform the path integration with 1-vertex 
boundary condition. 
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2.3. 1-Vertex Contribution 

We now calculate the contribution of a single 1-Vertex, assuming that this 1-vertex 
is at p point. Remember? In A field description, a single 1-vertex is a infinite small 
circle C around point p, with the inner of C cut out. And A field will increase 27r after 
counterclockwise walking along C a circle. This is the boundary condition of A path 
integration with a single 1-vertex located at p. One should integrate over p to calculate 
the total contribution of a single 1-vertex. 

We start from Z ( p.3|) , but now the path integral for A has the 1-vertex boundary 
condition, we denote this functional integral as Z[C]. 

Now will first integrate out A field in functional integration Z[C], all the calculation 
is the same as the previous considered case (without vextex), excerpt for the exponential 
factor exp f dB A dA\^ , now it will contribute a boundary term. 



exp 



^B{p) [ dA 
.^^ Jc 



(2.6) 



Since A will increase 27r around C, we can end the result 



Z'[C] = / [VB] exp[iB{p)] exp 



— / \\dB\ 



After integrating over p, we will get the contribution of a single 1-vertex. 

This calculation tells us that in B field description, vertex can be created out by using 
the local field operator exp[ii?(x)], this insertion of this operator will create a 1-vertex at 
X. Likewise, to create -1-vertex, one should use local operator exp[— ii?(x)]. 

We can now calculate out the total contributions of dilute gas of 1-vertex and -1- 
vertex, in terms of B field description. By noticing (^]^) we can see that the contribution 
of this dilute gas is a factor 



exp 



g I cos(S) 

R2 



(2.^ 



where the integration over comes from the integration over the position of vertexes 
(and anti- vertexes) . 

Now, in terms of B field description, the total partition function, include vertexes 
contributions (in dilute gas approximation), can be calculate as. 



ZI = / [DB] exp 



47r Jr2 



+ g I cos(S) 

'R2 



(2.9) 



where we use the symbol ZI to remind us that we have made 1-vertex and dilute gas 
approximation. 
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2.4- Quantum Turbulence 

We now turn to the BKT phase transition and quantum turbulence. From the ex- 
pression of Z/, one can see that: at extremely low temperatures, phase fluctuations are 
very large, they erased the contributions of vcrtexes. We can completely ignore the vertex 
at long distance. The spin correlation {J{x)J{0)) can be easily calculated with A field 
description 



Wile at high temperature, phase rigidity is large, vertex term contribution is very impor- 
tant, and B field are classically fixed at n. One can expanse B field around tt, to see that 
the system is gaped, the order is short range. 

The critical point can be determined by determine when the vertex fluctuation cos(S) 
becomes marginal. This can be seen by the standard fleld theory of B description 



When T > 1/4, cos(i?) is relevant, the system has energy gap, the renormalization 
group will flow to trivial at long distance. When T < 1/4, cos(i?) is irrelevant, the 
renormalization group will flow to a free boson CFT at long distance. When T = 1/4, the 
critical point, cos(i?) is a marginal items, when g = 0, the system will flow to a Z2 orbifold 
model of the c = 1 CFTs, and gcos{B) is the deformation along orbifold branch. 

As mentioned above, we have made 1- vertex and -1-vertex approximation. One should 
expect that di-vertex contribution will be a similar cos(2S) term. When the temperature 
is raised to 1, cos(25) will also become relevant. Therefore, our 1-vertex approximation 
only work when T is less than 1. When the temperature is high enough, all the vertexes, 
with different topological quantum numbers will become importance, di-vertex may be 
separated into two different 1-vertex, two may form a single di-vertex, etc., etc., such a 
process will continues and continues at long-distance. For example, if you stimulate a high 
topological quantum vertex. The vertex will separate and separate, all the way to, say, k- 
vertex, the fluctuations become relevant. At this point from 1-vertex to k- vertex, all their 
fluctuations, transformations, and interactions will be very important. Let us imagine, if 
k is a big number, large enough to make the k-vertex into a macro stuff, what happens? 
A little bit like the turbulence it? This is the quantum turbulence. 




T/2 



(2.10) 



1/(2T) = 2, T = 1/4. 



(2.11) 
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